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Scalar field dark energy evolving from a long radiation- or matter-dominated epoch has character- 
istic dynamics. While slow-roll approximations are invalid, a well defined field expansion captures 
the key aspects of the dark energy evolution during much of the matter-dominated epoch. Since 
this behavior is determined, it is not faithfully represented if priors for dynamical quantities are 
chosen at random. We demonstrate these features for both thawing and freezing fields, and for 
some modified gravity models, and unify several special cases in the literature. 
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I. INTRODUCTION 



The acceleration of the cosmic expansion can be ex- 
plained by a multitude of models suggested in the liter- 
ature but few or none derived from first principles. We 
look, therefore, for common characteristics among classes 
of models and for physically motivated generic behaviors. 
One example of effective characterization is the designa- 
tion of the thawing and freezing patterns of dark energy 
evolution [l[. These involve scalar fields leaving or ap- 
proaching the behavior of a cosmological constant in the 
dark energy equation of state w and its time variation 
w' — dw/dlna. Here we seek further unifying features in 
the evolution of the dark energy field. 

For inflation a slow-roll approach, neglecting higher 
time derivatives in the field evolution, can be used. This 
is not valid for dark energy, even at early times when w 
may be near — 1 and the dark energy is a small contribu- 
tion to the total energy density, and we must develop a 
different formalism. Several suggestions for specific ap- 
proximations exist in the literature and we will see that 
these can be unified into a single approach. This can 
also be extended in part to the freezing class of fields, 
traditionally difficult to characterize generically 

One key conclusion is that one must take into account 
that our universe is old: the scalar field has been evolving 
for many Hubble times in a background that was initially 
radiation-dominated then matter-dominated. This de- 
fines particular initial conditions and determines the dy- 
namical behavior. Employing random, Monte Carlo ini- 
tial conditions may lead to underrepresentation of thaw- 
ing or freezing behavior, due to neglecting the physics of 
a long evolution. 

In S|ll] we show why slow roll conditions are invalid 
for dark energy and discuss new methods for evaluat- 
ing the dynamics. We provide a specific example of a 
complete solution in mil\ Identifying a particular char- 
acteristic combination of parameters in ^IVl we show how 
the physics constrains the dark energy evolution within 
the classes of dark energy. We investigate extending this 
relation to some modified gravity scenarios in ifVl 



II. FIELD DYNAMICS 

Cosmic acceleration is given by the condition a > on 
the scale factor a, where dots represent time derivatives. 
For a state of near exponential expansion, the Hubble 
parameter H — d/a follows 



\H/H^\ « 1, 



(1) 



so the Hubble parameter is nearly constant. If the scalar 
field 4> provides the dominant contribution to the expan- 
sion, as it does for inflation, then the potential must 
be nearly constant, leading to an equivalent condition 
iV^/y I ^ 1 where V^ = dV/dcj) and we work in units with 
SttG = 1. This is often referred to as the (first) slow-roll 
parameter. Another implication is that the field acceler- 
ation (j) <C ZHcf) and so can be neglected in the equation 
of motion, or Klein-Gordon equation. 

However, these slow-roll conditions on T^ or rely on 
scalar field domination, and this is not valid for the dark 
energy field during its evolution, even today (since the 
matter density Vt,n is not negligible). In [3|, ratios of 
terms in the Klein- Gordon equation were defined as 



X 



H(j) 



Y = ^ 



-X 



Vi, X + 3' 



(2) 
(3) 



and shown today (when ri,„ k, 0.28) to be far from a 
slow-roll regime. In Fig. [T] we show explicitly that even 
for a thawing model that has equation of state today 
wq = —0.9, near the de Sitter w = —1 value, no term in 
the equation of motion is negligible compared to other 
terms, and also that \V^/V\ <C 1 does not hold, at any 
point in the evolution. 



A. General 

With all terms retained, the Klein-Gordon equation for 
the scalar field (p is given by 



3H^ = ~V^ 



(4) 



The ratios X, Y oi the terms are of interest for several 
reasons: 1) they indicate whether the potential driving 
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FIG. 1: Even a model evolving only from 'w{z ^ 1) = —1 to 
wo = -0.9 (here a PNGB model with / = 1, see qiIC|l cannot 
be described by a slow-roll formalism. Generically, none of the 
ratios of terms in the scalar field equation of motion, or the 
potential derivatives, is small enough to permit neglect of a 
term or to allow use of the slow-roll approximation. 



factor fli — > 0) 



) 

da' 



Q a' H{a') \ a 



- V^{a'). 



(7) 



Note that the ratio of the friction term and the driving 
term in the Klein-Gordon equation can be evaluated as 



3g0_ rda' fa'Y H{a)V^{a') 

V^ Jo a' \aj H{a') V^ia)' ^ ' 



It is also convenient to consider the ratio of the field 
acceleration to the driving term, 



Y= — 



-V^ ~ 3i/0 



/o a' \a 



3 r 



H{a) V^a') 
H{a') V^{a) 



- 1 



.(9) 



This is a formal solution because both H and V^ depend 
on (h itself. 



B. Asymptotic 



term or Hubble friction term is more influential (as used 
in []| to motivate thawing and freezing behaviors), 2) 
the ratios can be rephrased in terms of the tracks in the 
equation-of-state phase space w'-w, and 3) under certain 
conditions, such as tracking (constant w determined by 
the background expansion [3, 4]), X and Y will be con- 
stant. 

To begin, we obtain an implicit solution by mul- 
tiplying Eq. ((4|) through by an integrating factor 
exp{3 / dt' H{t')}, where the lower limit is arbitrary. In- 
tegrating from some early time ti to a later time i2 yields 



Some instructive cases can be evaluated directly. Con- 
sider a model where V^ is nearly constant in an epoch 
dominated by a component with equation of state Wb, so 
that H ^ Q-3(i-i-iu6)/2^ This is generic for the thawing 
class. Evaluating Eq. ([9]) gives 



.* oo — 



l + Wb 

3 + Wb' 



or Yoo = —1/3 for the matter dominated era, Y^^ 
for the radiation dominated era. 



(10) 

-2/5 



kt2) 



3f'^dt'H{t')l 



e -"2 






(5) 



Note that 



3 r*i dt' H(t') I I x3 

e ■''2 "• ' = [axlaiY , 



(6) 



so for t\ early enough relative to ti this factor is very 
small and the term involving 0(ii) is irrelevant. That is, 
the initial speed does not matter due to the high Hubble 
drag in the early universe. This leaves (taking the scale 



C. Thawing Evolution 

The evolution of the thawing field can be determined 
iteratively from Eq. ([7]) , taking into account that initially 
the matter contribution dominates in H and that the 
field moves relatively little from its initial value. Taylor 
expanding the potential around that initial field value 
(which without loss of generality can be set to zero), 



F(0) = ^, + K/,,0 + i^<^^,0' + 



(11) 



U' 



-\pmia) -f p^{a)\ = - {pra{a) + V; + . . .] (,12) 



where the subscript i denotes the value aX (j) = (j>i =0, pm 
is the matter density, and p^ is the dark energy density-'^ . 
At lowest order we find 



m 



2 V, 



(f)i 



[a = - 



"^ 



In second order 



2 K*,;a3/2 



y<pit 



(«) = -n 



9 Honlri^ 
2 V^, 



9pm(a) 



1-^ 

5 

-I 



2 V^, 
Q Pm{a) ' 




(14) 


2prn{a) 9pm{a)J 
6Vi J p,n{a)_ 


:i5) 

(16) 



Thus the criterion for the validity of the first order solu- 
tion for 4> is mainly that Vi/ pTn{a) <C 1, i.e. dark energy 
does not dominate, as expected. There are no slow roll 
V^/y ^ 1 or V,f„f,/V <C 1 conditions or other stringent 
condition on the potential derivatives. Indeed one could 
argue that having V and its derivatives be of the same 
order in Planck units is natural, as in the technically nat- 
ural PNGB model [| with V ^V^ [1 + cos(<?!)//)] where / 
is a symmetry energy scale of order unity. If V^^/F ^ 1, 
we see from Eqs. P3|) - (IT51) we might have to reevaluate 
the field expansion; see for example the next subsection. 
From these expressions we can calculate Y directly, 
substituting in the first order corrections to H , (j>, and 
V: 



Y{a) 



2 
15 



8 Vif,rj,. 

135 ~V^ 



V,, 



Pm{a) 



(17) 



The parameter Y starts off constant, at Kx,, only chang- 
ing as the scalar field rolls sufficiently far or its energy 
density starts to donunate the Hubble expansion. This 
behavior is evident in Fig. [TJ Also note that Y is not 
particularly small, and hence no one term dominates in 
the Klein- Gordon equation and ignoring <p is invalid. 

It is similarly straightforward to determine the lead- 
ing correction to the dark energy equation of state w by 
substituting the first order expressions into 



w 



Irk'^ 



<l?-v 



y^ + v 



to find 



VI- 



27 Vipmia) 



(18) 



(19) 



^ One can avoid specifying H explicitly by rewriting the Klein- 
Gordon equation as 



+ (1 + Q^ _ (1 + ^)n^) ^' + 3(1 _ n^) 



pm(a) 



0, 



where a prime denotes d/d\na. Since H^, 1 + w, and V^/ pm{<i) 
are all small at early times, this form is convenient for identifying 
the order of each term. 



The time variation of the equation of state is 

w' = 3(1 + 10), (20) 

the result obtained by [l[. That is, the evolution of the 
field in a matter-dominated universe fixes the asymptotic 
behavior of the dark energy for such thawing fields. 

D. Non-analytic Potentials 

As alluded to above, when the derivatives of the po- 
tential become large at a point, the Taylor expansion 
approach can break down. Consider a class of potentials 
with a singularity at some (/)^, e.g. 



V = 14-K(0-0.)" 



(21) 



with n non-integral and positive (we discuss negative n 
below) and cf) > 4>i,. This represents an inverted (concave) 
potential with the field rolling away from a maximum at 
(j) — (j)i, (eventually to negative infinity but that will not 
concern us regarding early time behavior) . 
We can find a thawing solution by trying 



Aa"" 



(22) 



in the Klein-Gordon equation during the matter- 
dominated era (easily generalized to other background 
evolution), or equivalcntly into Eq. ([7]). The result is 



A = 



Vn 2n(2 - n)2 
H^n^ 9(4 -n) 



1/(2-") 



(23) 
(24) 



The ratio of the field acceleration to potential slope 
terms, the equation of state, and derivative of the equa- 
tion of state become 



Y 

1 + w 
w' 



Hp^-m .2 2 37i/(2-n) 



3n 



(l + w) 



(25) 
(26) 
(27) 



as the field starts to roll. 

Note that for n > 2, in this Ansatz the field starts with 
a large kinetic energy, or equivalently w is positive and 
large, so we restrict to < n < 2 (for n = 2 the Ansatz 
fails). As n — > 2, the dark energy shoots away from 
w = — 1, acting more like sublimation than thawing. For 
the two integer values of n within this range, n — 0, 1, the 
potential has no singularities and these results agree with 
the Taylor expansion of the previous subsection. As well, 
if the field starts frozen away from the singularity then 
Taylor expanding the field works and the w' — 3{l + w) 
trajectory is the early-time solution. 



If instead we consider negative n, we can ignore the T4 
term at early times and (making the potential convex) 
we end up with a tracking field [O]- The equations for v, 
A, and Y above still hold but now 



Wn<0 



(28) 



and thus w' = 0. 



E. Unifying Relations 

Not only the dynamical trajectories but the relations 
between the dark energy density and the equation of 
state have characteristic behavior for each class of mod- 
els. Combining Eq. p9)) with the first order expression 
for the dark energy density, 



rirf, = 



V 



Praia) 



shows that 



n, 



271 



w 



A2 



(29) 



(30) 



where A — ~V^/V, here considered to lowest order in an 
expansion in 0. (Note that generally A^ = 2e, where e 
is the conventional first slow-roll parameter.) We return 
to this relation between the parameters in mVl now we 
simply explore some implications of the existence of such 
a relation. 

Note that to first order the relation between il^ and w 
does not depend on higher derivatives of V than the first 
derivative, and the relation between w' and w does not 
depend on V or its derivatives at all to this order. This is 
part of the unifying power of such an analysis, that any 
scalar field dominated by the background Hubble friction 
must behave in a simple, determined manner. 

We can now compare our result to other approaches in 
the literature that assumed specific potentials or approx- 
imations. The simplest case is the limit V = constant or 
A ^ 0. This is of course just the cosmological constant 
and the equation of state never leaves w — ~1. (Skat- 
ing models [3, Q have a large initial motion ^i but this 
quickly redshifts away, as a~^, and the field comes to rest 
at w — —1.) Next is the linear potential [9|, [lO|, where 
V^ — constant, discussed in the next section. In general, 
though, potentials will have higher order derivatives that 
are not zero or depend nontrivially on the first derivative. 

Thawing models have been studied with approxima- 
tions, such as taking A — constant (tur ning the exponen- 
tial potential's tracking behavior [ll|, \1^ into thawing 
by starting it from a frozen state) but approximating fJ^ 
or w 13]. Indeed [l3| noted a version of the relation 
((30)) then holds asymptotically. Another parameter in- 
vestigated in a first order expansion about a constant 
value is k = — A/(l + X/3) [l3|. Explicitly incorporating 
Eq. (pOl) . [15] expanded in the energy density fl^ about 



the asymptotic solution to form an "algebraic thawing" 
model, which is actually valid to second order. Interest- 
ingly, current data show the algebraic thawing model is 
statistically a better fit than ACDM [16^ . 

All of these cases follow the unifying first order solu- 
tion ([30]) but diverge at higher order. Each one basically 
chooses different ways to truncate or close the hierarchy 
of higher order equations. To understand how quickly de- 
viation from the unified solution, or of the validity of the 
field expansion approach, occurs, we note it requires par- 
ticular combinations of Q^ and 1 -I- w to be much smaller 
than one. So we do not expect these analytic field solu- 
tions to be valid up to the present. Nevertheless, in many 
cases they are good approximations until surprisingly re- 
cent times; e.g. for the model of Fig. [1] the relationship 
([501 holds to 3% (8%) until z = 2 (1). We discuss this 
further in i ^IVI 



III. DARK ENERGY DENSITY AND THE 
LINEAR POTENTIAL 

In addition to understanding the dark energy equa- 
tion of state, we are often interested in the observables 
directly, such as the Hubble parameter or the distance- 
redshift relation. These involve the dark energy density, 
given by the sum of the potential V{4>) and kinetic energy 
0^/2 discussed in the previous section. 

The dark energy density at some epoch relative to its 
current value is given by the ratio 



? 



,{a) ^ '^^' + V 
'0.0 ^'J'I + Vq 



(31) 



where the subscript indicates current values. The cur- 
rent value for the equation of state is 



Wo 



\4>l - ^0 



(32) 



Differentiating the dark energy density with respect to 
time and using the Klein-Gordon equation, we find 



di _ 3iJ(/)2 1 - Wo 
di ~ Vq 2 



or 



d£, _ 3(/)M - Wo 
da Vo 2a 

To remove the explicit appearance of </> we use 

3i/0 



1 + Y 



Y4> ' 



and to eliminate Vo we employ 

<^g^ ^0.o(l + yo)^ 

Vo mlvo 



2(l + u;o) 



Wo 



(33) 



(34) 



(35) 



(36) 



with the result 



da 






1 



Wo i?o 



y,' 



H' VI, 



1 + Y{a) 

l + Yn 



(37) 



Note that the Hubble parameter is given by 

i/V-ffo = f^mo a-^ + (1 - r!,„o) ^(a) (38) 

so the density ^(a) is an implicit function only of the scale 
factor a - except for the dependence on the potential 
slope V^{(j){a)). This dependence occurs for Y{a) as well 
through Eq. ^. 

In the special case of a potential linear in (/>, the quan- 
tity V^ is constant. This was treated by Q as one of the 
first dark energy models and more recently as a textbook 
example by 17[ • Then the two equations p7|) and ([9]) in- 
volve only the independent variable a and the dependent 
variable ^. (This will be a good approximation for any 
potential with a slope varying with a sufficiently slowly.) 
Despite being coupled integro-differential equations, they 
are actually simple to solve numerically. One starts with 
an initial approximation for the Hubble parameter with 
£, set to unity for all values of a, that is initially H{a) is 
appropriate to a cosmological constant. In this case, two 
parameters fix the solution: Oa = 1 — ^mo and wq. 

The value of Y{a) is not known until the solution is 
found iteratively, but for the first iteration the cosmolog- 
ical constant Ansatz 



-0.1 




FIG. 2: The acceleration response function Y{a) for any thaw- 
ing field has common characteristics. Here we plot the case 
of a linear potential with wq = —0.777, fi^^o = 0.76 (compare 
to the PNGB case shown in Fig. [TJ . The solid curve is the 
exact solution, while the dotted curve shows the analytical 
approximation for wo = —1, O^^o = 0.76, which is used as the 
starting point for the iterative solution. 



H{a)/Ho = ^{l-nA)a-^ + nA = h{a) 



(39) 



gives 



1 + Y(a) 



1 — r^A h{a) 



h{a) 



+ - 



2f] 



1/2 



In 



/i(a)2 
h{a) 



flA 



n 



1/2 



h{a) + Q 



1/2 
A J 



(40) 



The dark energy density calculated with the first order 
quantities is quite accurate, even when wq is not very 
close to — 1 , and convergence upon iteration is extremely 
rapid. By taking a further derivative of ^, one can also 
derive the relation for the time variation of the dark en- 
ergy equation of state followed by the linear potential, 
w'q = dw/d\na{a = 1) « —1.3(1 -I- wq). Results found 
iteratively for the linear potential are shown for Y(a) and 
the density ^(a) = Pcj,{a)/p4,.Q in Figs. [5] and [31 



IV. FLOW PARAMETER 



As we saw in J jIIEl evolution away from the frozen 
state involves a synchronized deviation in both Q^ and 
l + w away from zero in a particular proportional or scal- 
ing manner that persists to quite late times. This im- 
poses a constraint that phenomenological models must 
obey, and parameter priors or Monte Carlo scans must 



take into account. Other forms of dark energy, such as 
modified gravity, might also follow similar scaling rela- 
tions (see, e.g., [l^) and we investigate this in the next 
section. Such constraints are important as we lack first 
principles theories for dark energy, as scalar fields or mod- 
ified gravity. 

Motivated by the relation in Eq. ([5(I)l . we define a new 
quantity we call the flow parameter. 



F : 



1 



nA^ 



(41) 



where A = — V^/y, fi^, and w are all functions of time. 
This is related to the friction vs. potential slope term 
ratio 1 -I- y by 

1 



F^-il-wYil + Yy, (42) 

and to the phase space evolution of dark energy as 



y' = -3(l-w2)[l-(3F)-i/2] 



(43) 



These relations follow generally from the Klein-Gordon 
equation and from the definition of Y, for all a and with- 
out any requirement for matter domination. We note 
that F holds close to its high-redshift, asymptotic value 
for a substantial part of the evolution, even when fiw, is 
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FIG. 3: The dark energy density relative to the present, ^(a), 
changes httle for thawing fields. The curve shows the solution 
for the linear potential with wq = —0.777, fl^,o = 0.76; note 
the density used in the first iteration (dotted curve in Fig. ^ 
is just the cosmological constant behavior, unity for all scale 
factors. By contrast, for a typical tracking field ^{a ^ 1) 2> 1, 
e.g. ^{a = 10"^) « 10* for V ~ (^"^ 



FIG. 4: A particular combination of dark energy parameters, 
F = {1 + w)/[i}4,{dV/d(l}f], has fixed high redshift behavior 
for each class of models, 4/27 for thawing and 1/3 for freezing 
(red, dotted lines). Also, it is nearly conserved during the long 
evolution until dark energy begins to dominate. 



not much less than one. Figure [4] shows the evolution 
of F for thawing and freezing (in fact modified gravity) 
examples. 

As seen in Eq. (j30p . for any thawing model evolving in 
a matter dominated universe the field flow begins with 



F = 4/27 (thawer) . 



(44) 



This is equivalent to the condition w' = 3(1 + z«) at early 
times, but is preserved for more of the evolution. The 
other major class of scalar field models is tracking mod- 
els. These have constant equation of state at early times 
in a background dominated universe, and from Eq. 
this implies 



F = 1/3 (tracker). 



(45) 



This can be viewed as Y{a <C 1) = (1 + u'oo)/(l — Woo), 
where Woc is the high redshift dark energy equation of 
state. The non-analytic thawing solutions with < n < 1 
interpolate between the regular thawing case (cf. n = 1) 
and tracking models {n < 0): 



F 



4 /2 
3 



4- n 



(46) 



Again, there is a substantial range of validity for F be- 
ing constant even as fi^ grows to an appreciable fraction. 



This relation imposes a particular high redshift behav- 
ior on a whole class of phenomenological models. Ef- 
fectively, these flow behaviors define physical priors that 
must be included in any Monte Carlo simulation of scalar 
field evolution. As shown, they hold until quite recent 
epochs, z ^ 1 — 2, not just for a ^ 1. Assuming instead 
a random dynamical state, rather than that determined 
during the many matter-dominated e-folds of expansion, 
is equivalent to a fine tuning of the scalar field initial con- 
ditions to avoid the natural evolutionary path. See [l^ 
for an interesting analysis of the effect of the difference in 
priors between those uniform in w'-w and those uniform 
in a version of the field parameters. 



V. MODIFIED GRAVITY AND EXPANSION 

The accelerating expansion of the universe could indi- 
cate a deviation from general relativity rather than the 
presence of a new scalar field. It is of great interest to 
explore beyond scalar fields and see whether the evolu- 
tion of expansion and matter density growth parameters 
probes modified gravity. For example, [18, 19] motivated 
a deviation from the general relativistic growth behavior 
scaling as r2^(a), or a^ in many cases, while [20| assumed 
a variation as a, and [2l| left it as a fit parameter a". 
As seen in [2l|, the observational consequences for the 
deviation index s are significant. This is too large and 



complex a subject to address here, but we examine how 
some modified gravity theories respond to the parame- 
ters treated above for scalar field explanations of dark 
energy. 

For gravity extended beyond Einstein relativity, the 
parameters $7^, w, V, etc. are effective quantities but 
we can still define them in terms of modifications of the 
Friedmann equation from the general relativistic, matter- 
dominated case. 



A. DGP Braneworld and H" 

We first consider the Dvali- Turner [22| modification 
(also see l2l|) 



H' 



8TrGpmia)/3 + {I - r!„,o) H^ {H/HoT, (47) 



where flm,o — ^m(a = !)■ This was motivated by con- 
sideration of extradimensional theories, with the index 
a depending on boundary conditions between our four 
dimensional universe and a higher dimensional bulk vol- 
ume. The DGP braneworld cosmology [23, l2a] is the 
special case a — I. Following \2^ , this acts as an effective 
freezing scalar field with 



w 



2-a 



-n.. 



w' — 3w(l + w) 



l--(l + u;) 
a 



(48) 
(49) 



This provides the opportunity to follow the flow parame- 
ter for a freezing (indeed tracking) field starting far from 
w = —I, rather than for a thawing field. From Eq. (|^5)) . 
we find the effective flow parameter 



F(a) 




f72 



2{l + a-^)^n^ + {^y 



r(i50) 



Considering the early time limit of F, we expect in 
the matter-dominated universe that Qm —^ Ij so F —i- 
1/3. That is, the flow parameter indeed agrees with the 
tracking value of Eq. (|^5)l . Asymptotically, the effective 
equation of state is constant at w = — 1 -I- a/2 and one 
can in fact show that the effective potential has the form 
of an inverse power law potential V ~ </)"" with index 
n = 2a/{2-a). 

In the future, the H'^ term will come to dominate (as- 
suming a < 2) so that Vt,n ^0. Its effective equation 
of state then approaches w = —1. The flow parameter 
F -^ 4/3, which corresponds to w' — 'iw{\ -t- w). This 
flow parameter is the maximum value for freezing scalar 
flelds and the minimum value for barotropic fields (see 
[2g | for the barotropic case). In this limit, the i?" effec- 
tive potential approaches a quadratic form, V ~ T4o [1 + 



(3/8) (0OO — 0)^], where the field has an asymptotic max- 
imum value ^oo and V^ = 3H^ [1 - 0„(z = 0)]2/(2-"). 
Note that this class of models thus effectively incorpo- 
rates a cosmological constant. Figure \5\ shows the full 
solution for the effective potential. 
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FIG. 5: The effective potential of the modified gravity case 
involving H" is plotted in units of Hq for a = 1/2, 1. At high 
redshift (small <j)) it possesses the flow parameter of a tracking 
scalar field and acts like an inverse power law potential, while 
in the far future it freezes to a cosmological constant state, 
approaching it as a quadratic potential. 



B. f{R) Gravity 

Another class of theories where dark energy is an effec- 
tive quantity arising from extending gravity is f{R) the- 
ories, where the action involves a function of the Ricci 
scalar, here considered in addition to the usual linear 
term (so f — corresponds to general relativity). [231 
gives the modified Friedmann expansion equation and 
an expression for the effective dark energy density pde- 
We consider the deviation from the early, high curvature, 
matter-dominated era and write the equation in terms of 
the effective total equation of state of universe wtot , 



/" 



1 3 

2 + 2'"*°' 



W]f 



- -(1 - 2w;tot - iw^ot + w'tot) / 

= 2A7rGpd,il-2wtot~3wl^ + wU, (51) 

where W = {3-3wtot-'^5w'^ot-^Wtot+^wiot+^wtotw[at- 
w[',,)/{l-2wtot-3wl,+wi^,). 
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Now we consider the evolution of the departure from 
general relativistic matter domination with total effective 
equation of state Wtot — 0. Since the total equation of 
state Wtot — wVlwio-), where the dimensionless effective 
dark energy density ^w{a) = (87rG/3)pdo/^^, then we 
have flw{a), wtot, and its derivatives all of the same order. 
The dark energy density fi^u (a) <C 1 at the epoch consid- 
ered so we can expand in this small quantity. The zeroth 
order solution, with wtot = = Qw{a), gives f ^ a^ 
with p = (-7 + V73)/4, as found by ^. This homo- 
geneous solution acts as an initial condition that quickly 
becomes unimportant. The first order departure from 
matter domination and general relativity has 



fi{a)=AH^{a)ny,{a), 



(52) 



where A is of order unity, and shows that gravity deviates 
from general relativity at the same rate as the effective 
dark energy density evolves (cf. [la])- The proportional- 
ity quantity A does not have to be constant, in general, 
just of order unity (thus any time variation is on Hubble 
scales or longer). 

If we restrict A to be constant, then we can solve the 
equation in terms of effective dark energy equation of 
state and phrase this as the flow parameter 



Fia) 



3{l-w^)^ 



[3A- 



(5/2)H 



(53) 



We can find a tracking solution with w' = 0, yielding 
A = — 3/[l — {5/2)w— 3w^], so a given choice of constant 
A corresponds to a given constant w. Since the dark 
energy density evolves as pdc ~ a~'^'-'^+™) and the Ricci 
scalar is dominated by matter pm ''^ a~^, in this limit 
the gravitational modification looks like f{R) ~ i?^+™. 
However, considerable freedom exists to choose other so- 
lutions, e.g. with A varying. 

Note that the physics governing the true scalar field 
evolution can be very different from that operating for 
modified gravity, so there is no expectation that the same 
relations should hold. Flows unlike the well-determined 
quintessence behavior may provide hints of modified 
gravity. 



VI. CONCLUSIONS 

Our universe is old, having expanded by a factor of per- 
haps 10^* since the last period of acceleration during the 



inflationary epoch. This strongly affects the evolution of 
a scalar field that may give rise to the current epoch of 
acceleration and determines some key properties of the 
dark energy. Although conventional slow roll approxi- 
mations are invalid, we show that analysis in terms of an 
integral relation between the Hubble friction and poten- 
tial driving terms and a well characterized field expansion 
can give insights into the evolutionary behavior. 

For the case of thawing fields, the field expansion pro- 
vides a clear initial track in the w-w' phase space, a uni- 
fication of a number of interesting special cases, and a 
rapid convergence in the evolution of the dark energy 
density. For tracking fields, the ratio between the Klein- 
Gordon terms reaches a constant value. In both cases, 
the evolution of the deviation of Q^ from zero and of 
the tilt 1 + w scale in a manner constrained by the long 
matter (and radiation) dominated era. Phrased in terms 
of a flow parameter combining the scalings, this ratio is 
nearly conserved until quite recent times, z « 1 — 2, when 
the dark energy flnally begins to take over. 

This physical behavior means that dark energy dynam- 
ics is not random, or equally probable, e.g. in the sense 
of a uniform prior over w-w' , but is focused - "flows" - 
in specific ways. We have also tested this for two modi- 
fied gravity theories and found some similar behavior but 
also some deviations that could offer clues to the nature 
of the acceleration. While one can always arrange ini- 
tial conditions such that the dark energy comes out of 
the matter-dominated era with arbitrary behavior, this 
involves fine tuning. The oldness of our universe does 
provide a natural path for dark energy dynamics. 
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